Splitting in the Fermi surface of ZrTe3: a surface charge density wave system 
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The electronic band structure and Fermi surface of ZrTes was precisely determined by linearly 
polarized angle-resolved photoelectron spectroscopy. Several bands and a large part of the Fermi 
surface are found to be split by 100-200 meV into two parallel dispersions. Band structure cal- 
culations reveal that the splitting is due to a change of crystal structure near the surface. The 
agreement between calculation and experiment is enhanced by including the spin-orbit potential in 
the calculations, but the spin-orbit energy does not lead to a splitting of the bands. The dispersion 
of the highly nested small electron pocket that gives rise to the charge density wave is traceable 
even in the low-temperature gapped state, thus implying that the finite correlation length of the 
long-wavelength modulation leads to a smearing of the band back-folding. 



I. INTRODUCTION 

In a one-dimensional metal the atoms form chains and 
specific electronic states are confined to travel predomi- 
nantly along these chains. This leads to a Fermi surface 
(FS) topography that is characterized by parallel sheets 
as the band dispersion is significant only along the chain 
direction. The corresponding FS nesting makes the sys- 
tem susceptible to instabilities such as spin or charge 
density waves (SDW and CDW), if assisted by a coupling 
to the respective magnetic or lattice degrees of freedom. 
Chain-like arrangements are found in a number of crystal 
structures and CDW phenomena are observed in varied 
systems such as blue bronze, KCP, (TaSe4)2l and the 
prismatic chains of NbSea and TaSs.^ The ideal CDW 
should be an isolated Peierls chain, but in a crystal a cer- 
tain coupling between parallel chains is always present. 
Reduced coupling of adjacent chains is achieved in surface 
CDW systems such as In/Si(lll)2 and Au/Ge(001).3 

To date it is still very difficult to predict whether a 
material will exhibit a charge density wave. Fermi sur- 
face nesting is certainly an important ingredient and can 
be studied by ground state electronic structure calcula- 
tions of the unmodulated structure. The formation of a 
statically modulated CDW ground state requires a cer- 
tain coupling between the chains so that below the Peierls 
transition a three-dimensionally ordered ground state can 
exist. Phonons eventually destroy the CDW order at the 
Peierls transition temperatures, but on the other hand 
the electron-phonon coupling EPC is also needed to as- 
sist the formation of a static distortion where the lattice 
degrees of freedom follow the electronic modulation. It is 
thus highly interesting to study how the CDW reacts to 
slight changes in the chain structure and coupling, which 
can be induced by pressure,!^ strain, or, in the case dis- 
cussed in this paper, by the presence of the surface. 

We present an extended data set of photoelectron spec- 
troscopy data of the layered material ZrTes. A cleaved 



surface of a van der Waals bonded layered system is usu- 
ally assumed to have a bulk-like termination and the 
absence of dangling bonds makes the formation of sur- 
face states unlikely. Thus the surface-sensitive photo- 
electron signal is often regarded as representative of the 
bulk electronic structure. Indeed, in ZrTes a previous 
photoelectron spectroscopy stud}'^ found the momentum 
structure of the charge density wave and its temperature 
dependence in good agreement with the complementary 




FIG. 1: (a) Projection of the unit cell onto the a-c plane of the 
structure showing the chains in the Te layer adjacent to the 
van der Waals gap. (b) First Brillouin zone of the monoclinic 
structure and surface Brillouin zone (SBZ). (c) Schematic rep- 
resentation of the two main sheets of the Fermi surface with 
the nesting vector. qcDW represents the projection of qcDW 
onto the SBZ. 
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and less surface sensitive optical conductivity data.'^ We 
find, however, that our data are not compatible with the 
expected features of the bulk electronic structure and 
conclude that the signal, as well as the results of other 
studies,^ ^ ^ represent a special surface system that is 
modified from the bulk structure, with the interesting 
consequence that this modified system shows potential 
differences to the related bulk of ZrTea. 

In the bulk, the CDW transition of ZrTea is well char- 
acterized. It was first observed as a resistivity anomaly 
at TcDW = 63 K for measurements along the crystal 
a and c-directions.^ The periodic lattice distortion was 
observed by diffractioiP^ and by a giant-Kohn-anomaly 
behaviour of a soft phononf^ Below Tqdw the system 
is left in a metallic state, where superconducting fila- 
ments emerge at very low T (Tc = 2 K).^^ The open- 
ing of a gap in the electronic structure does not affect 
the entire Fermi surface. In fact, the Fermi surface 
consists of at leas t two principal sheets (Fig. [ijc)) as 
seen in calculationJ^^ and experiments.^ ^ A central 
3D Fermi surface is derived from holes that propagate 
mostly along the prismatic (ZrTeajoo chains along the h- 
direction. The elongated cross section in the a*-6* plane 
allows for a certain nesting along 6* . A much better nest- 
ing condition is found in the small electron pocket of the 
Te bpx band which propagates along the Te-Te chains 
formed by the lateral bonding of adjacent prisms. This 
quasi one-dimensional (qlD) Fermi surface with principal 
propagation direction along a was identified to match the 
observed modulation vector qcDW = (0.07, 0, 0.333).^^ 
Strictly speaking, the qlD FS consists of two sheets due 
to the inequivalence of the two Te atoms in the Te-Te 
chains. In the out-of-plane direction c the structure is 
much more weakly bonded since only van der Waals 
forces act between layers. The band dispersion is con- 
sequently much weaker. Some dispersion has, however, 
been predicted by all calculations. In particular the qlD 
Fermi surface can be seen to follow the inclination of the 
Brillouin zone boundary (BZB) given by the monoclinic 
P2i/m structure.S^ 

The Te layer adjacent to the van der Waals gap is 
of special interest. It is formed by the inequivalent but 
almost equidistant Te*^^^ and Te*^^^ atoms and its band 
structure has been discussed in Ref. 12, Due to the weak 
bonding it is the most likely surface termination layer 
in a cleaved sample. The nesting qcDW of the artificial 
two-dimensional single layer would be along the F-B di- 
rection. 

This paper is organized as follows: in the next sec- 
tion we detail the experimental procedure and compu- 
tational details. Sect. IIIll shows the measured data and 
discusses the observed band dispersions and splittings. 



Sect. |IV| discusses the implications of the observed split- 
ting through comparison with calculations. Sect. [V| re- 
ports on the observations concerning the CDW gap. The 
conclusions are summarized in the final section. 



II. EXPERIMENT AND CALCULATIONS 

Single crystals of ZrTea were grown by vapour-phase 
transport in the form of large shiny platelets. Photo- 
emission data were acquired from samples freshly cleaved 
at room temperature by peeling off the top layers with 
Scotch^^-tape. In an ultrahigh vacuum of typical resid- 
ual pressure 2 • 10~^^ mbar the samples remained clean 
and sharp spectra could be acquired over a two days pe- 
riod following cleavage. The sample normal and high 
symmetry directions were determined from the shape of 
the sample, the direction of the fibres along the b axis and 
the symmetry of photoemission angular scans. The Fermi 
level and the total energy resolution {AE = 17 meV) 
were determined from the metallic Fermi cut-off of a poly- 
crystalline gold film. 

The angle-resolved photoelectron spectroscopy 
(ARPES) experiments were performed at the linearly- 
polarized undulator beamline BL-1 at the Hiroshima 
Synchrotron Radiation Centre.^^ The experimental 
station is equipped with a low-temperature two-axis 
goniometer. The angle-dispersive detection direction 
of the electron spectrometer is in the plane of linear 
polarization of the photons. When a mirror plane of 
the sample is aligned with this plane, the spectra will 
only contain contributions from orbitals that are even 
with respect to the mirror plane (selection rule). When 
the sample is tilted about a horizontal rotation axis, 
both even and odd states may be detected depending 
on the magnitude of the tilt. Fermi surface maps have 
been acquired by sequentially tilting the sample across 
and away from this symmetric situation to cover the 
irreducible portions of the two-dimensional in-plane 
Brillouin zone. Additional data were acquired using a 
similar set-up with even higher angular resolution at 
the Swiss Light Source using elliptically polarized light 
(Figs. [4] and 0). 

Mapping of the photoelectron angular distribution to 
momentum k\\ was performed according to standard par- 
allel projections. The locus of the Brillouin zone bound- 
aries was determined from the symmetry of the data sets. 
While ARPES could in principle detect the full three- 
dimensional momentum vector j -*^^ ' -*^^ ' we find that the ob- 
served band dispersions agree perfectly for data sets ac- 
quired at different photon energies hu = 36.7 — 58.8 eV. 
Variations are seen only in the relative intensity of fea- 
tures, but not in their binding energies at corresponding 
Our data therefore display a two-dimensional band 
structure and we denote the high-symmetry points in the 
surface Brillouin zone as r, Y, B and D (cf Fig. [ijb)). 

Two distinct sample mountings were employed that 
lead to different orientations of the polarization vector 
with respect to the sample. With the 6-axis mounted 
horizontally the central 3D Fermi surface was accurately 
mapped since the tilt corresponds to the dispersion along 
F-B. No mirror plane of bulk ZrTea is present in these 
scans, but states that are odd in a are suppressed in this 
sample orientation. With the 6-axis mounted vertically. 



3 



the qlD Fermi surface and the B-point (where the two 
principal sheets approach) are accurately mapped. Here 
the selection rule implies that at the B point only even 
states with respect to the a-c plane are observable. 

Peak positions, intensities and widths were extracted 
from the measured data wherever stable peak fitting 
could be performed using model functions consisting of 
two or four Lorentzians. In the case of momentum dis- 
tribution curves (MDC) at constant binding energy, the 
widths of the Lorentzians were restricted to be equal 
for equivalent bands and a constant background was in- 
cluded in the fit function. Typical widths were of the or- 
der of 0.07 A (0.03 A for the higher resolution data 
sets). In the case of energy distribution curves (EDC) at 
constant position in momentum space, the resulting sum 
of Lorentzians and a linear background was multiplied 
by a Fermi distribution function with a pre-determined 
Fermi energy and an effective temperature (width of the 
Fermi cut off) that takes into account the finite energy 
resolution. Typical widths of EDC features ranged from 
70 meV up to 120 meV for very steep bands, where 
the momentum resolution has an effect on the spectral 
widths. 

Ab initio band structure calculations were performed 
using the Wien2k code.^^ For bulk ZrTea the crys- 
tal structure from the literatur^^ was optimized until 
the forces on all atoms converged below a limit of 2 
mRy/Bohr. The augmented plane wave basis was ex- 
panded to a wave-number limit kmax given by r • k^ax = 
7, where r = 2.5 Bohr is the radius of the atomic spheres. 
The total energy was integrated over a momentum space 
mesh ofl3xl9x7 points which was reduced according 
to the crystal symmetries. For a calculation in a slab ge- 
ometry the unit cell was doubled in the c-direction and 
a vacuum layer of 10 Bohr thickness was inserted. The 
two-dimensional electronic structure was calculated on a 
26 X 40 X 1 mesh and the atomic positions were relaxed to 
less than 2 mRy/Bohr resulting in slightly expanded lay- 
ers. Other calculation parameters were identical to the 
bulk calculation. Fermi surface cuts were interpolated 
from the Fermi level crossings of all bands in a regular 
40 X 20 grid in the symmetry-reduced part of the first 
Brillouin zone. 



III. BAND DISPERSIONS AND SPLITTINGS 

Dispersion maps at two photon energies across F in the 
a-c mirror plane are shown in Fig. [2] The photoemiss- 
ion intensity is plotted on a linear grey scale versus the 
binding energy Eb and the wave vector along F-B. 
The two data sets differ in sample temperature T and 
photon energy hu. Further data sets have been acquired 
in the photon energy range hu = 36.7 — 56.8 eV. As 
the total momentum of the photoelectrons is varied with 
their energy a signature of the dispersion along /c^^ ^ is 
expected. However, except for variations in the relative 
intensity, the observed bands, their binding energies and 
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FIG. 2: ARPES dispersion maps along the direction F — B 
measured with photon energies (a) hu = 42.5 eV and (b) 
hi/ = 58.8 eV at low temperatures as indicated. Dots in the 
left section of (b) are peak positions from peak fitting. Dashed 
lines mark the observed dispersions in selected regions as a 
guide for the eye. The top of each panel shows the intensity 
distribution at Ef with the peak positions marked by arrows. 
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FIG. 3: ARPES dispersion map along the direction F — Y 
measured with hu — 42.5 eV and T = 30 K. The top section 
shows the intensity distribution at Ef- Dots are peak posi- 
tions from fitting as shown in the right section for an EDC at 

r. 
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dispersion match between all the data sets. Thus the 
data represent the two-dimensional electronic structure 
of a ZrTes layer. 

The top section of each panel shows a cut along kx 
(MDC) at Ep. The strongest feature is a small pocket 
^^a^' at the B point just below Ep- A second Fermi surface 
feature ^^bi shows two Fermi crossings at kx = 0.32 A 
and 0.4 A , together with its symmetry-related copies. 
It consists of two parallel dispersing bands with a split- 
ting of ^ 100 meV. The third feature "c" is a completely 
occupied band with a band apex at F. Two copies with 
parallel dispersions may again be distinguished, although 
the relative intensities differ strongly. The left section of 
Fig. shows an EDC at F, where the two peaks due 
to this band are highlighted by the resulting Lorentzian 
functions from peak fitting. The peak positions from 
this fitting procedure are marked as dots in the panel. 
The splitting is ~ 100 meV at F and increases to almost 
200 meV at the lower band apex. 

|0A-1 I0.24A-1 I0.8A-1 
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FIG. 4: Fermi surface map of ZrTes as measured by ARPES 
at /iz/ = 45.1 eV and T = 200 K. The global scale is given and 
the crystal symmetry is indicated by the 2D Brillouin zone 
boundaries. The features marked "a", "bi" and "62" corre- 
spond to Figs. [2] and [3] The two panels have been measured 
with linearly polarized light in different directions as indi- 
cated by arrows. Markers show the position of MDC peak 
fits. Horizontal lines (a-f) correspond to cuts shown in Fig.|5] 
and arrows at the top to Fig. [7| The inset shows a higher res- 
olution data set around B acquired using elliptically polarized 
light. 

An orthogonal map through F in the direction of Y is 
shown in Fig. [S] Feature "c" is again visible as a double 
peak (EDC in the right section of the figure). The most 



prominent feature ^^62 is a steeply dispersing band with 
a Fermi crossing at /c^ = 0.15 A . Due to the strong 
intensity differences of the two parallel dispersing copies, 
its splitting is revealed only by careful peak fitting as 
shown on the top part of the figure for an MDC ai Ep 
and by tracing the same band along the Fermi surface. 

The Fermi surface map (Fig. |4| was measured by ac- 
quiring dispersion panels at various tilt angles away from 
the high-symmetry plane. The intensity at Ef was ex- 
tracted by integrating over a window of 40 meV around 
Ep. The two panels of Fig. [4] were thus measured with 
orthogonal photon polarizations as indicated. The data 
were not symmetrized or distorted. At T = 200 K, the 
electronic structure is still unaffected by the CDW tran- 
sition^ ^ and the Fermi surface is representative of the 
unmodulated high-temperature unit cell. 
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FIG. 5: (a-d) Momentum distribution curves at Ef for cuts 
parallel to ky at various kx as indicated by horizontal lines in 
Fig. [4] Solid lines are the resulting curves from peak fitting, 
(e) and (f) ARPES dispersion maps parallel to ky at kx as 
indicated. The peak positions from peak fitting are marked 
by dots. Map (e) coincides with curve (d). 

The locus of features ^'a", ^'bi^' and ^'b2^' is indicated 
in the Fermi surface map in Fig. [4j The identification 
of these features enables us to attribute the small pocket 
^'a^' to the small electron pocket of the qlD Fermi surface. 
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The intensity is seen as a wide band on both sides of the 
BZB. Features "bi" and "b2" form the more corrugated 
3D Fermi surface. 

The pocket of intensity of the qlD sheet (feature "a^') 
is visible ah along the line B - D. The centre between 
the two MDC peak positions forms a slightly curved line 
although it should follow the straight line of the BZB. 
We assign this artificial curvature to slight motions of 
the beam spot as the sample is rotated. The observation 
of a dispersion (see also Fig. [t] below) allows the deter- 
mination of the two peak positions of the Fermi wave 
vectors from the MDCs. The inset displaying high res- 
olution data around B reveals that at least one Fermi 
surface sheet in the qlD pocket does not reach all the 
way to the B point as a closing arc is visible. However, 
a high intensity is seen in a region everywhere along the 
BZB including B. Inside the sharp Fermi arc, towards 
the edges of the plot inset, a further Fermi surface sheet 
is observed which yields intensity in a sharp line following 
the BZB. 

The 3D sheet ( ''hi and ^^62 '0 has the shape of a fiat- 
tened oval. The splitting leads to two concentric copies. 
The Fermi wave vectors have been determined by peak 
fitting four Lorentzian peaks to the MDCs at Ep and 
the resulting peak positions are indicated by markers in 
Fig [4j The momentum splitting at location ^^62 is re- 
duced to almost zero. This is due to the steep disper- 
sion of the band as seen in Fig. [3) The splitting be- 
comes readily visible in its evolution along the Fermi 
surface as shown in Fig. |5] The MDCs shown in the 
top sections (a-d) are well described by peak fitting four 
Lorentzian peaks of equal width. In panels (e) and (f) 
showing maps along ky for values of kx where the FS 
sheet closes, peak fitting was performed on the EDCs, 
identifying two parallel copies of the band and allowing a 
splitting of 120 — 180 meV to be determined. Note that 
panel (f) corresponds to the occupied part of feature "bi 
in Fig. [2j The intensities vary due to the different direc- 
tions of the polarization of the light. In particular the 
contrast to another band obscuring the view in Fig. [2] is 
greatly enhanced. 



IV. ORIGIN OF THE SPLITTING 

The splitting of both the band at F ^'c^^ and of the 
3D Fermi surface by 100-200 meV is not expected in the 
bulk electronic structure. Both bands are derived from 
highly directional Te p-orbitals with some hybridization 
with Zr d states. The low-symmetry crystal structure 
does not permit a degeneracy of these bands (except for 
the Kramers degeneracy due to the inversion symmetry, 
which forces spin- up and spin-down to be degenerate). 
Possible explanations of a splitting of this order of mag- 
nitude include (z) ferromagnetic symmetry breaking lead- 
ing to an exchange splitting or (ii) spin-orbit (SO) split- 
ting of the kind found in the closed p-shells of atoms. 
The former {i) would require the sample to be ferromag- 



netic, which has not been observed. The latter (ii) was 
proposed by Pacile et al. in the related ZrSes, HfSes 
and ZrSs.l^ In these semiconducting systems a band at 
r is observed in ARPES spectra that disperses like fea- 
ture ^^c^' in our data and shows a similar splitting. The 
arguments in favour of SO splitting are (a) the order of 
magnitude of the splitting which matches atomic expec- 
tation for the chalcogenide p-shell and (b) an increase of 
the splitting in the selenide compared to the sulphide. 
This splitting would require, however, that the underly- 
ing band has a degeneracy which could be lifted by the 
SO-interaction. Other origins of the splitting could be 
of a geometric nature, such as {Hi) bi-layer splitting or 
(iv) inhomogeneity of the surface leading to a superpo- 
sition of features from differently cleaved or clean and 
dirty portions. We exclude this last possibility (iv) since 
the splitting and intensity ratios of the features were ob- 
served reproducibly on various cleaves in different sample 
orientations and in different laboratories. 

In order to further identify the character of the ob- 
served bands, we performed a first-principles band struc- 
ture calculation. The result of the Fermi surface and 
band dispersion of bulk ZrTea in the a*-6* plane are 
shown in Fig. [6ja) and (c). The agreement of this cal- 
culation and the experiment is rather poor, although we 
may assign most of the bands. Feature "d^' at F and 
the 3D Fermi surface are readily identified as the band 
at 0.5 eV binding energy (in agreement with the data) 
and the large hole pocket extending in the F-B direc- 
tion, respectively. The qlD Fermi surface is found over 
a large portion along the BZB from B-D and consists 
of two sheets due to the two inequivalent atoms which 
form the Te-Te chains. The band ^'c^' is assigned to the 
hole pocket around F, thus incorrectly crossing the Fermi 
level. Such a small pocket with few carriers could easily 
be an artifact of the calculation and has little significance 
for the total energy. The calculations were performed us- 
ing a scalar-relativist ic treatment of the SO interaction, 
thus any degeneracy-lifting due to the SO energy should 
be visible in the calculation. Overall, the agreement of 
this theory with the experiment is so poor that even a 
precise assignment of all bands is hard to achieve. 

An inspection of the calculations shows that the bands 
in question are non-degenerate (except for the Kramers 
degeneracy). The splitting due to SO interaction {ii) 
must therefore be ruled out. Of the scenarios listed 
above we are left with {Hi) a bi-layer splitting due to 
"almost but not quite" equivalent atoms in the effective 
unit cell. A surface reconstruction which could lead to 
an effectively larger surface unit cell and hence inequiv- 
alent atoms has not, however, been observed by LEED. 
We propose to consider an out-of-plane relaxation lead- 
ing to an effective bi-layer situation near the surface. To 
model the surface we construct a repeated slab of two in- 
equivalent ZrTea layers. The slabs are separated by a 10 
Bohr thick vacuum spacer and the structure is relaxed in 
Wien2k. This leads to a widening of the layers and the 
two layers move slightly apart as shown in the schematic 
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FIG. 6: Fermi surface and band structure calculations of (a-b) bulk ZrTes and (c-d) a two- layer slab of slightly distorted 
ZrTes. The unit cells for the ab-initio are shown schematically. In the slab unit cell, the displacements with respect to the bulk 
configuration are indicated by arrows. 



unit cell in Fig. [6[b). The electronic structure in this 
stable configuration is strictly two-dimensional. The real 
surface configuration of a cleaved sample is probably close 
to the bulk configuration up until the second-outermost 
layer. The outermost layer, deprived of its neighbours, 
can relax slightly into the vacuum and expand. The 
probing depth of ARPES at Ek = 45 eV is likely to 
show only features from the outermost or two outermost 
layers. This surface configuration could in principle be 
modelled by a large slab calculation comprising many un- 
reconstructed and reconstructed layers. However, such a 
calculation is beyond the scope of this paper. We believe 
that we capture the most important features with our 
simple bi-layer model. 

The calculated Fermi surface (FS) and band structure 
of the bi-layer are shown in Fig. |6jc) and (d). The 3D 
FS sheet is split into two concentric sheets. The qlD FS 
sheets now consist of four sheets, all of which are parallel 
to B-D and fulfill good nesting conditions. The small 
hole pocket around F still appears in the calculation, but 
the dispersion reveals that only one of the split bands is 
slightly unoccupied. Except for this mismatch of binding 
energy, the overall agreement of this calculation with the 
experiment is much better. In particular, feature ^'d^' and 
the hook-like feature ^^e^' close to B are recognized in the 
data (cf. Fig. |2|. The splitting of the 3D Fermi surface 
is well-reproduced in this calculation. The splitting of 
feature ^'bi at F is seen in the calculation, although it is 
smaller (50 meV) than that observed experimentally (100 
meV). In agreement with the experiment, the splitting 
increases to about twice this value along the dispersion 
F-B. 



Based on the agreement of the two-layer slab calcula- 
tion with our experimental data, we consider the surface 
of cleaved ZrTes to be strongly relaxed in the out-of-plane 
direction. The detailed atomic structure cannot be de- 
rived from our spectroscopic data, but the observed split- 
ting suggests a configuration with a reconstructed surface 
layer which is almost equivalent to the bulk layers. This 
conclusion leads to the interesting consequence that other 
photoemission studies on these surfaces (refs. 5.7 8) may 
also have probed a special surface configuration of the 
material. 



V. THE CDW GAP 

Peierls-type instabilities are released by a momentum 
dependent gap in the electronic structure ai Ep corre- 
sponding to the incommensurately modulated structure. 
The gap of ZrTes has previously been observed,^ ^ but 
only as a reduction of spectral weight at the Fermi level 
within the intensity on the BZB. The dispersion of the 
electron band could not be resolved. The higher resolu- 
tion of the spectrometers used in the current study allows 
the band to be seen as a parabolic dispersion along kx for 
all sections along B-D. As the temperature is lowered 
and spectral weight is shifted away from the Fermi level, 
the band dispersions remain traceable up to although 
the intensity is reduced within the gap. Fig.^shows three 
sections along kx at B (equivalent to Fig. |^a)), near B 

at ky = 0.24 A and at D. At least one electron-like 
band is observed at B, two at /c^ = 0.24 A ^ and again 
one at D. 
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FIG. 7: ARPES dispersion maps at three different position along the Brillouin zone boundray BZB from B to D at low T = 30 K 
(e-g) and T — 200 K (h). The Fermi wave vector on either side of the BZB is determined by peak fitting as shown in (a-d) of 
MDCs at Ef- In panels (g-h) the dispersion is traced by fits to the MDCs by two Lorentzians of equal width. The intensities 
of these fits are shown to the right of the panel together with a Fermi occupation function (dashed line) at the corresponding 
temperature. The half-height of the intensity curve at T = 30 K (h) corresponds to the magnitude of the gap Sg and is marked 
by an arrow. The data sets (a-b) and (e-f) were acquired with a higher resolution spectrometer and elliptically polarized light, 
whereas (c-d) and (g-h) were acquired with linearly polarized light. 



The Fermi surface is thus partially gapped, leaving am- 
ple sheets for both metallic conductivity below Tcdw 
and superconductivity. The 3D sheet does not show any 
change with temperature, the qlD sheet has a fully de- 
veloped gap in a large region around D. On approaching 

B, the gap is reduced and no gap is seen for ky < 0.2 A . 
The inset of Fig. [41 shows a closing arc of one sheet of the 

o — 1 

FSat/c^ = 0.16A . AtBa high intensity is, however, 
still present. We assign this intensity to another FS sheet 
that is not identified by the calculations. This FS sheet 
at B deserves special attention as Yokoya et al.^ have 
found a spectacular increase in spectral weight close to 
Ep at low temperature in this position. In our data the 
MDC at Ep of the intensity at B is compatible with a 
single Fermi surface crossine: (Fig.jTfa)). The band is not 
observed when the photon polarization is along 6* (see 
Fig. [4|, thus it is even with respect to the 6*-c*-plane. 
It was suggested that a van Hove singularity (vHs) could 
be present at this point Ep due to the f act that the 
3D and qlD FS sheets touch each other.'^^'^ Any involve- 
ment of the 3D FS sheets can be ruled out, however, since 
its Fermi wave vectors are found to be /c^^ = 0.32 A and 
kx = 0.4 A respectively for the two split sheets, at a 
safe distance from the B-point (kx = 0.55 A ). Our cal- 
culations fail to predict a Fermi surface sheet at B and 
cannot assign the observed band. The hybridization of 
the Te p^^-band with the band dispersion from feature 
^^c^' could still lead to a small vHs (cf. Fig. [6]^d)). 

In the region ky > 0.2 A ^, where the CDW gap 
was observed in previous studies,!^ the band dispersion 
and intensity distribution is determined by a peak fitting 
analysis of the MDCs. In a large region around D, the gap 
is fully developed at T = 30 K.f^ Surprisingly, the band 
dispersion remains unambiguously traceable up to E^i? as 
shown in Fig. [7|^f-g) . The shift of spectral weight away 



from Ep is seen as a linear increase of the band intensity 
over a binding energy interval of 150 meV. This compares 
to a constant intensity over the whole binding energy 
range derived from a similar analysis at room tempera- 
ture (Fig. [7|^h)). The expected back-folding of the band 
due to the new periodicity is not observed and the gap 
value is determined as the binding energy where the band 
acquires half of its maximum intensity Eg = 65± 10 meV. 
An upper band apex (whose binding energy would cor- 
respond to the gap value) is not found. 

ZrTea is a rather special CDW system since the com- 
ponent qcDw''^*-' = 0-07 along a* is very small and corre- 
spondingly the wavelength of the modulation along the 
chains is very long (A 80 A). This property is con- 
served in the cleaved surface of ZrTes and we determine 
2k p = 0.06 A = 0.11- a*, slightly larger than sug- 
gested by the bulk superstructure. This nesting suggests 
a CDW in the surface layer as (g^,^^) = (0.11,0). The 
long wavelength and the finite correlation length of the 
Peierls modulation, even well below the transition tem- 
perature Tcdw = 63 K, lead to a smearing of the band 
folding in momentum space. Since the Peierls distortion 
only creates a slight modulation in the underlying crystal 
structure, it is not surprising that the back-folded band is 
not visible in the data. However, the removal of spectral 
weight from the Fermi level is still assigned to a partial 
CDW gapping of the Fermi surface. 

VI. CONCLUSIONS 

We have presented a detailed Fermi surface map and 
band dispersion study of ZrTes. The basic features of 
the Fermi surface, i.e. a large 3D body around F and 
a qlD sheet along the Brillouin zone boundary B-D are 
consistent with previous results, but the 3D Fermi surface 
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as well as a band at F are found to be split by 100-200 
meV. This splitting is explained by a surface relaxation 
leading to a bi-layer splitting of the electronic structure. 

The qlD FS sheets are seen all along the BZB B-D. 
The band is traceable by peak fitting with a Fermi wave 
vector 2kF = 0.06 A . This opening of the electron 
pocket is seen everywhere along the BZB, although at 
least one sheet of the Fermi surface closes in an arc 
at /c^ = 0.15 A . A second sheet is identified inside 
the nested sheet and at D they are degenerate. Very 
good nesting conditions are therefore found all along B- 
D. The CDW gap is identified as a shift of spectral 
weight away from Ep m. gap oi £g = 65 ± 10 meV 

for ky > 0.2 A . This gap is not seen as a back- folded 
band, due to the long wavelength and finite correlation 
length of the Peierls distortion. 

This leads to the conclusion that the surface of ZrTea, 
as well as other MX3 systems, does not have a simple bulk 
termination. Even though the crystal probably cleaves in 
the van der Waals bonded gap between Te layers and the 
surface termination has no dangling bonds, the outer- 
most layer relaxes into a configuration which is slightly 
different from the situation in the bulk. The electronic 
levels of this surface system split in a fashion similar to a 
bi-layer splitting. All observations by highly surface sen- 
sitive techinques, in particular the VUV ARPES studies 
performed by Yokoya et al.^ Starkovicz et al. and Facile 
et al. (Refs. 5 , 7 and 8), and probably also the studies on 
NbSea by Schafer et al^ have therefore probed a spe- 
cial surface system of reconstructed MX3 chains, which 



closely resembles the bulk but is slightly modified. Since 
a small modification of the atomic arrangement can dras- 
tically infiuence the CDW, this opens up the possibility 
to deliberately manipulate the CDW material on its sur- 
face and thus study the CDW phase transition in new 
and controlled systems. 
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